published in: Phys. Rev. Lett. 101, 260503 (2008) 



Bell test for the free motion of material particles 
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We present a scheme to establish non-classical correlations in the motion of two macroscopically 
separated massive particles without resorting to entanglement in their internal degrees of freedom. 
It is based on the dissociation of a diatomic molecule with two temporally separated Feshbach 
pulses generating a motional state of two counter-propagating atoms that is capable of violating 
a Bell inequality by means of correlated single particle interferometry. We evaluate the influence 
of dispersion on the Bell correlation, showing it to be important but manageable in a proposed 
experimental setup. The latter employs a molecular BEC of fermionic Lithium atoms, uses laser- 
guided atom interferometry, and seems to be within the reach of present-day technology. 
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Introduction. — The possibility of entangling macro- 
scopically separate, non-interacting particles challenges 
our classical view of the world by putting into question 
the concepts of realism and locality [IJ. Nowadays en- 
tangled states are routinely established with photons y . 
Also the internal states of material particles have been 
entangled, e.g., using non-classical light as a carrier for 
the quantum correlations [3|, 4,|5L|6|, |7|, or their Coulomb 
interaction in an ion trap a ISl- However, the origi- 
nal discussion of entanglement focused on the motional 
state of massive particles, whose spatial separation is a 
dynamic feature of the entangled two-particle wave func- 
tion 10|,lll|. The latter is spatially extended, unlike with 
internal entanglement, where the positions only play the 
passive role of separating the parties. Since the positions 
and momenta are quantum observables with a direct clas- 
sical analog, an observation of non-classical correlations 
arising from macroscopically distinct phase space regions 
would therefore be a striking demonstration of the failure 
of classical mechanics. 

A convincing demonstration of non-classical correla- 
tions between two parties requires the experimental vi- 
olation of a Bell- type inequality [ij]. In the simplest 
case, it involves detecting dichotomic properties on each 
side, such as the polarization of spin-1/2 systems, where 
a maximal violation is obtained if the spins are in a Bell 
state, say |$) = (| T)i| T)2 + e'"^! i)i| i)2)/\/2. Such a 
dichotomic property is not readily available in the mo- 
tional state of free, structureless particles. One possi- 
bility is to consider observables which have no classical 
analogue, such as the displaced parity or pseudospin op- 
erators used to discuss entanglement in the 'original EPR 
state' 13[ |l4|. While they are expedient for character- 
izing the quadratures of light fields, their experimental 
implementation seems exceedingly difficult in the case of 
free, macroscopically separated material particles, where 
only position measurements are easily realized. 

Here we aim instead at producing a motional analogue 
of the Bell state |$) and implementing a Bell test which 
requires a simple position measurement in the end. The 
idea is to expose a diatomic molecule to a sequence of two 
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FIG. 1: (a) The two particles in a DTE state are each charac- 
terized by an early and a late wave packet component result- 
ing from two different dissociation times, (b) Each particle 
is passed through an unbalanced Mach-Zehnder interferom- 
eter with a switch deflecting the early wave component into 
the long arm, while conducting the late component through 
the short arm. (c) The path difference is chosen such that 
the early and the late wave packets overlap. Detecting the 
particle in an output port, at a given phase ip and splitting 
ratio -&, amounts to a dichotomic measurement analogous to a 
spin 1/2 detection in an arbitrary orientation. The existence 
of DTE is established, implying a macroscopically delocalized 
two-particle state, if the correlations at the single-particle in- 
terferometers violate a Bell inequality. 



temporally separated dissociation pulses. Each of the two 
counter-propagating dissociated atoms then has an early 
and a late wave packet component corresponding to the 
two possible dissociation times, and their correlation may 
be called "dissociation time entanglement" (DTE), see 
Fig. [T] (a) . If one regards these components as spin state 
analogues in the motion (say, early corresponds to spin 
up, late to spin down), the DTE state shows the same 
structure as the Bell state |<I>). Switched, unbalanced 
Mach-Zehnder interferometers on each side then serve 
to mimic arbitrary spin rotations, and the detection of 
each particle in one of the interferometer output ports 
completes the Bell test, see Fig. [T](b,c). 

The DTE state is a variant of the "energy-time entan- 



gled" state introduced in [l5|, and it is closely related 
to the "time-bin entanglement" of photons [2, uM, Il7| . 
which has been used e.g. for establishing non-local cor- 
relations over fiber distances of more than 50 km 'l8] in 
a similar setup as in Fig. [1] As a main difference, the 
DTE state is not composed of two single-particle prod- 
uct states, but it superposes their relative coordinate, 
while the center-of-mass state remains separable during 
the two-pulse dissociation process. 

Our use of DTE reflects the necessity to come up with 
a state generation scheme appropriate to material parti- 
cles. Their finite mass and internal structure entail sub- 
stantial complications which require a careful investiga- 
tion of whether a Bell violation can be expected at all. 
Most prominently, one must account for the unavoidable 
wave packet dispersion, and only the recent progress in 
manipulating ultra-cold molecules (such as their conden- 
sation [19| and controlled dissociation [20|) suggests the 
possibility to generate motional states that allow one to 
keep the detrimental effect of dispersion under control. 
We note that other ways to demonstrate nonlocal cor- 
relations of molecular dissociation products have been 
proposed in 2lU22|. 

We will show below that, by appropriately choosing a 
sequence of magnetic field pulses, the Feshbach-induced 
dissociation of an ultracold ^Li2 molecule within a guid- 
ing laser beam can generate a motional state of macro- 
scopically entangled atoms that is capable of violating 
a Bell inequality. But first, to clarify the implications 
of dispersion, we determine the correlation function for 
a generic DTE state subject to correlated single-particle 
interferometry. 

The general DTE Bell test. — It is natural to take the 
bound molecular two-particle state to be separable in 
the center-of-mass (cm.) and the relative (rel) coordi- 
nate. Denoting by r the period between the two dis- 
sociation pulses and assuming the transverse motion to 
be frozen in the ground state of the guiding laser beam, 
the longitudinal part of a pure DTE state then takes 
the form |^dte) = (ui°l|*o) + c'^-|*o))/V2, where 
l^o) = |V;§-'°-)(|V'5°') -hP|Vo''))/V2. It involves the free 
time-evolution operator U^J, the parity operator P, and 
a relative phase 4>t determined by the details of the two- 
pulse dissociation process. The center-of-mass state of 
the original molecule |'0o'") is taken to rest in a wave 
guide, while I'^Jj'^') propagates into positive direction. 

Indicating the output ports of the two interferome- 
ters J = 1,2 by (Tj- = ±1, the immediate experimen- 
tal outcome is characterized by the two-time probabil- 
ity density pr(cri,CT2;ii,i2) for detecting particles be- 
hind the respective ports at the times ti,i2- It de- 
pends on the phase settings ipj and mirror transmis- 
sion angles dj of the interferometers, and the disper- 
sive evolution of the DTE state leads to a complicated 
fringe pattern as a function of ii, ^2, and t. How- 
ever, as suggested by the above analogy with a dis- 



crete Bell test, a robust quantity characterizing entangle- 
ment is obtained by integrating the port-specific proba- 
bilities over all times, P(ti.,(T2 = /pr(o'i, 172; ii,i2)dtidt2. 
For any reasonable model of the time-of-arrival detec- 
tion, the correlation function can be equally evaluated 
by means of the projections IIo-i <8 Ho-a = 11^2 onto the 
(unbounded) spatial regions behind the respective ports, 

P^,,^, = limt^oo ||n^jU^,t|^DTE)|P. 

Since we are only interested in the interferometer out- 
put, the effect of dispersion is best incorporated by time- 
dependent scattering theory, which separates the 'raw ac- 
tion' of the interferometers, described by the S-matrices 
Sj , from the free dispersive time evolution U^ I . The pro- 
jection of the DTE state component |'0o'™')l''/'o°') onto a 
particular output-port combination ai , CT2 then takes the 
form 
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Here we assume t to be sufficiently large, so that the 
wave packets have passed the interferometers, and we 
take the entrance switches to be in the "on" configura- 
tion, i.e., routing towards the long arms. Implementing 
the phase shifts by the arm- length variations ij , the pro- 



jected S-matrices Si.™ — Ilcr^Sj 



are given by S^. ^^^^ 



eyip{ipj£j/ti) cosdj and S)j°.'^_i = eyip{ipj£j/h)sm'dj. 
For the "off" configuration we have correspondingly 
Sg,°J^]^ = siwdj and Sg.°J_i — —cos-dj. Analogous rela- 
tions hold for |V'§'")P|V'"'). 

The setup requires the dispersion-induced broadening 
of the wave packets to remain much smaller than the 
separation between the early and the late components, 
so that the switches can be changed in between. In this 
case the correlation function Pax,iJ2 can be evaluated by 
using Eq. ([T]) and its variants even for non-pure and non- 
separable initial two-particle states, |^o)(^o| — > Qa- One 
obtains 
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2)'^, and pr(pi,p2) = 



with p = ipi,P2) , i = {ii 
{pi,P2\go\pi,P2) the momentum distribution function. 
For simplicity we take here the beam splitters to be sym- 
metric {i)j = 7r/4) and the particles to have equal mass 
m. Note that ([2]) is independent of the total time of 
flight and invariant under momentum phase transforma- 
tions (pi,p2|*o) -^ exp[i^(pi,p2)](pi,P2|*o), which in- 
cludes spatial translations, thus rendering the correlation 
^^1,(72(^1,^2) a robust signal. 

For generic Gaussian states in the center-of-mass and 
relative motion Eq. ([2]) can be evaluated in closed form. 



The variances Cp ^^ 
teristic dispersion times, T^ 



and (Tp j.pj then determine charac- 
2TO?i/(Tp (.^ and T^ci = 



m?i/2(Tpj.^[, indicating the transition to a dispersion- 
dominated spatial extension of the wave packets. The 
expectation value of the relative momentum po.roi = 



mUfci/S defines the reduced wave length Arei = S/po,rei, 
which sets the scale for the non-local interference fringes 
in the explicit correlation function, 
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These DTE correlations can violate a Bell inequal- 
ity. This is seen from the structural similarity of ([3]) to 
the correlation function of the standard spin-l/2-based 

setup, Plf^^{ipi,ip2) = {l + criCr2C0s((pi -(^2)}/4, 

where the CTj = ±1 denote the spin measurement out- 
comes for analyzers tilted by the angles ipj with respect 
to the Bell state quantization axis. It follows from this 
analogy that an unambiguous demonstration of entangle- 
ment requires the fringe visibility of the correlation signal 
to exceed l/\/2 over at least a few periods. 

The dispersive suppression of this fringe visibility is de- 
scribed by those terms in ([3]) which depend on the charac- 
teristic times Tc.m. and T^ci- Specifically, the dispersion- 
induced distortion between the early and the late wave 
packet components is reflected in the Lorentzian reduc- 
tion factor and in the quadratic compression of the fringe 
pattern, while their envelope mismatch causes the Gaus- 
sian suppression. Based on this result one finds that non- 
local correlations can be observed, for Aroi/(''"Uioi) <C 1, 
provided (l + T2/r<2„, )(l-hT2/T^2j) < 4. In the following, 
we present a conceivable scenario for the generation of a 
DTE state, which meets these conditions. 

Experimental scenario for a DTE Bell test. — We sug- 
gest to use a dilute molecular Bose-Einstein condensate 
(BEG) produced from a 50:50 spin mixture of fermionic 
^Li as a starting point. It can be prepared efhcicntly 
and with near-perfect purity, displaying a huge lifetime of 
more than 10 s due to the Pauli blocking of detrimental 3- 
body collisions [23j . A truly macroscopic time separation 
between the two pulses, say r = 1 s, is thus conceivable, 
and for a realistic dissociation velocity of Uici =1 cm/s the 
de Broglie wave length Aioi = 13.3/im is large enough to 
pose viable stability requirements for the interferometers. 

The BEG is prepared in a red-detuned, far off-resonant 
laser beam (transverse trap frequency loq/2'it — 300 Hz), 
strong enough to guide the dissociated atoms towards 
the single particle interferometers. The intersection with 
a second, weak laser beam creates an elongated dipole 
trap for the BEG within the laser guide [24|. At the 
end of the preparation steps it is very shallow (depth 



number of molecules (on the order of 10 ) remains in 
the BEG. These can be taken to be non-interacting, so 
that the initial longitudinal center of mass state JV-'t) of 
the molecules is straightforwardly defined by the trap 
parameters. 

Each interferometer (with path length difference 
TUrci/2 ~ 5 mm) is implemented by two more red- 
detuned laser beams crossing the guide in a triangu- 
lar arrangement at small angles. While the crossings 
act as beam splitters, the required atom mirror may 
be realized using an evanescent light field or a blue- 
detuned laser beam perpendicular to the interferometer 
plane [25|, |26|]. The time controlled appliance of such 
perpendicular blocking beams could also implement the 
switch. However, a simplified setup could replace the 
switch by an ordinary beam splitter, at the cost of 50% 
post-selection. 

The fluorescence detection of the slow, strongly con- 
fined atoms at the two output guides can be done 
with near unit efficiency and with single particle reso- 
lution, since no particular spatial or temporal accuracy 
is needed. This single-particle resolution is crucial since 
events where more than one molecule gets dissociated 
are disregarded in the present scenario, relying on post- 
selection. In a more refined setup, it is conceivable to 
use a specially prepared optical lattice where each site is 
occupied by at most one molecule [27]. 

All this implies that the molecular dissociation in 
presence of the wave guide must meet a number of 
rather restrictive criteria to render a demonstration of 
macroscopic entanglement possible. Only the transverse 
ground state of the guide may be populated to admit the 
above quasi-one-dimensional description of the interfer- 
ometers, while both the momentum spread of the wave 
packets and the dissociation probability must be suffi- 
ciently small. In order to judge the feasibility of our 
experimental scenario, we now discuss how the dissoci- 
ated part of the state, |$bg(0): depends on the magnetic 
field pulse sequence and the resonance parameters. 

A Green function analysis within the two-channel 
single-resonance approach shows that after an arbitrary 



magnetic field pulse sequence (close to an isolated reso- 
nance) |<&bg(0) is described, for low energies, at positions 
far from the dissociation center, and at large times, by 

the asymptotic form |$bg(t)) ^ a^\v-of-)H%)V^°l\'^z), 
where U^, l is the free propagator in the longitudinal 
direction and where the transverse motion is frozen in 
the harmonic ground state, |(/9§™) and Ifl^})), respec- 
tively, of the guiding laser beam. The longitudinal 
state is determined by {pc.m.,Prci\^z) = C(p^^ /Amh + 
Prci/™^+2wG)(Pc.m.|V'T)/||C'||, whcrc C'(w) is the Fourfer 
transform of the closed channel probability amplitude 

C{t) and \\C\\^ = /dpe.m.dprcl|C'b?.m./4™?i + p2^l/TO?l-H 

2ujQ,)\^\{pc.m.\'<pT)\^ ■ C{t), in turn, is determined by the 
coupled channel dynamics as induced by the externally 
controlled magnetic field B{t). The dissociation proba- 
bility is given by |CbgP = WGabgA*rosAi?rcs||C'|p/7r?i2. It 
involves the background scattering length abg, the reso- 
nance width Ai?ios, and /Xics, the difference between the 
magnetic moments of the resonance state and the open 
channel. 

While the association of the molecules is best done 
at a broad resonance [23], shifting to a narrow reso- 
nance (e.g., ABics = ImG, ^rcs = 0.01 /xb, Obg = 
lOOao) one can ensure by choosing field pulses B{t) with 
short duration (e.g. T — 60 ms) that about a sin- 
gle molecule dissociates on average. A sequence of two 
square pulses with base value Bo and height Ai? sweeping 
over the resonance position B^^s then generates a DTE 
longitudinal wave packet of the required form, \^ z) = 
[Ui°| -h e''^-]|*o)/\/2. The momentum distribution is 
given by |(pc.m.,Prci|*o>P = PoP^sinc^ [(p^^ /4 
p2)/Ap2]|(p,.„.|VT)PA[b?.m./4 + P?, 
where we define Po/m = ^res{Bo + AB — Bres) — 2C/t — 
hiOG, as well as p'^/m = /XresAS and Ap^ = 2mh/T 
[28| . Instead of evaluating Eq. ^ directly with this mo- 
mentum distribution, it is more transparent to approx- 
imate the latter by Gaussians, which allows one to ap- 
ply the analysis following Eq. ([3]). For reasonable pulses 
Bq + AB — i?ros = 350 mG these Gaussians are centered 
at po.roi = ±?7iWi.ci/2 and po.c.m. = 0, with spreads of 
Cp.rei — l-l96mh/(jioT) and (Tp,c.m. = ^hw^mj^^ re- 
spectively. This yields Trei = 3.4 s and T^.-m. — 0.64 s, 
implying a visibility of about 72% in the correlation sig- 
nal, which exceeds the threshold value of l/\/2- 

Our analysis thus shows that the observation of a 
macroscopic DTE state is feasible with material par- 
ticles, even though dispersion poses tight constraints. 
The technological challenge is substantial, but not insur- 
mountable. Stable lasers are required and the magnetic 
pulse sequence must be reproducible with a relative ac- 
curacy of 10~^ from shot to shot, so that the relative 
phase between the early and late components, given by 

^r — [^UtT - fircsABT + ^res(Sros " Bo)T]/h + UJgT, 

varies less than 50mrad. Realistic choices of the laser 
wave length (1 ^m) and the vacuum pressure (10~* mbar) 
suffice to suppress decoherence due to scattering of off- 



cm./ ^ ' Prel 
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resonant photons or background gas particles. 

A great advantage of this setup is that no interferomet- 
ric stability is required between the two interferometers, 
so that truly macroscopic separations are feasible. More- 
over, the DTE state reveals its entanglement robustly 
since neither a spatial nor temporal resolution is required 
in the detection. 

This work was supported by the DFG Emmy Noether 
program. 
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